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ABSTRACT 

The present cosmic rate of Type la supernovae (SNela) suggests that about 
6% of all stars in binary systems with primaries in the initial mass range 
3 — 9 M© end up as SNela. If that is confirmed, the unavoidable conclusion is 
that SNela can only be explained by the single degenerate scenario. At most 1% 
of stars in binary systems in the above range end up as CO + CO WD pairs, 
with total mass equal to or larger than the Chandrasekhar mass. Given that the 
number of mergers from pairs of CO + He WDs that reach the Chandrasekhar 
mass is even lower, the conclusion strongly favors binaries containing just one CO 
WD as the progenitors of SNela, since the SNela production efficiency (relative 
to the instantaneous star formation rate) predicted for double degenerate (DD) 
pairs lies more than 3o" below the observational data, and the DD scenario can 
be rejected at more than 99% confidence level. Only if the SFR measurements 
from z ~ 0.1 to 2; ~ 0.5 are being underestimated by a factor of 6 while SNela 
rates are not, can we escape the above conclusion. We evaluate the numbers 
and characteristics of double WD systems with different chemical compositions 
(CO and He WDs) that should form and compare them with the observations, 
in order to check our predictions. Our conclusions appear robust after that test. 

Subject headings: cosmology: supernovae: general — binaries: close binaries — 
stars: white dwarfs 
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1. Introduction 

The evolution of intermediate and low-mass stars in binary systems is a complex 
subject that has, however, enormous importance in astrophysics. Of their way to the final 
end, a bulk of information is available from observations on the processes undergone by 
binaries as they evolve, i.e. common-envelope phases, accretion, outburst. 

A relevant new piece of information comes from the consideration of the number of 
Type la supernova explosions (SNela) that occur per unit mass spent in forming stars in 
different environments. So far, the exploration of the issue was restricted to nearby galaxies, 
where the combination of SNela statistics with star formation rate (SFR) history is poorly 
known. 

The approach explored here is to take the average of star formation history in redshift 
space and use the information on cosmic SNela rate to infer how many stars in binaries 
from low to intermediate mass end up as SNela. This indirect approach started a few 
years ago (Ruiz-Lapuente, Burkert, & Canal 1995; Ruiz-Lapuente, Canal, & Burkert 1997; 
see also Madau, Delia Valle, & Panagia 1998; Sadat et al. 1998; Ruiz-Lapuente & Canal 
1998; Yungelson & Livio 1998, 2000; Kobayashi et al. 1998; Dahlen & Fransson 1999; 
Ruiz-Lapuente, Casse, & Vangioni-Flam 2000) and it can help us to determine the nature 
of the progenitors of SNela. 

We evaluate the number of SNela exploding per unit comoving volume, in redshift 
space, and relative to the mass going into forming stars in the Universe. Despite 
the uncertainties still involved, the last years have nevertheless brought some crucial 
information. 

As we will see, if the efficiency of binaries in ending as SNela is 6% of the stars between 
3 — 9 M0, no other alternative than the single degenerate scenario appears as a reasonable 
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candidate to explain those explosions. The modeling on which that conclusion is based 
stands comparison with the statistics of the white dwarf (WD) pairs actually observed. 

2. The star formation rate and the SNela production efficiency 

In the last few years we have started to learn about the cosmic history of the star 
formation process, and several groups (Madau et al. 1996; Steidel et al. 1998; Hughes et 
al. 1998; Madau, Pozzetti, & Dickinson 1998; Blain et al. 1999) have begun to derive the 
global star formation rate p^{z). At the same time, the first measurements of the global 
SNela rates, ^ia{z), are being performed (Pain et al. 1996; Hardin et al. 2000; Hamilton 
1999; Hamuy & Pinto 1999; Pain et al. 2000). The SNela rate is doubly related to the 
SFR: through the time delay between formation of the SNela progenitor systems and 
explosion, and through the fraction of stars (binaries in this case) that give rise to SNela. 
The comparison between p*{z) and ^ia{z) thus contains key information on the nature of 
the so far elusive SNela progenitors. 

The work of Madau et al. (1996), based on the UV luminosities of galaxies in the 
Hubble Deep Field, complemented that of Lilly et al. (1995, 1996) at optical wavelengths 
in showing that the star formation activity steadily increases with z from the local, 
present-day Universe, up to 2; ~ 1.5. The original claim that p*(-2) peaks there to fall again 
at higher z has later been revised, in particular on basis to data at long wavelengths that 
would indicate that the star formation regions at high z are enshrouded in dust (see Hughes 
et al. 1998, for instance). p*{z) would level-off beyond z ~ 1.5 — 2, up to much higher 
redshifts. Observations in the submillimeter range (Blain et al. 1999), ISO measurements 
of the extragalactic background light (Rowan-Robinson et al. 1997; Flores et al. 1998; 
Elbaz et al. 1999), and the FIRAS and DIRBE experiments on board of COBE (Dwek 
et al. 1998; Fixsen et al. 1998), do in principle lift the veil on the star formation activity 
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in dust-obscured regions, but there is no redshift identification of tlie emission detected. 
Measurement of tlie SN rates at high z then would help to trace the star formation history. 
SNela are the most luminous ones. Thus, knowledge of the characteristic time delays 
between formation of their progenitor systems and explosion is most relevant. 

The measurements of the cosmic evolution of the SNela rate now extend up to 2; ~ 0.55, 
and there are already preliminary results for z ~ 1.1. Current results are summarized 
in Table 1. By combining our knowledge of K/a(-2) with that of p^{z), we can now gain 
new insight on the nature of the SNela progenitor systems. We introduce for that the 
"efficiency" £sNeia{z) of SNela production referred to the SFR, both per unit of co moving 
volume: 

SsNeia{z) = ^ia{z) yr'^ Mpc-^ / p,{z) Mq yr~^ Mpc-^ (1) 

The quantity £sNeia{z) (given in the fourth column of Table 1) is thus just the number 
of SNela per unit mass spent in forming stars at a given z, and it is independent from 
the cosmological model assumed. This efficiency reflects the evolutionary time scale (from 
birth to explosion) of the progenitor systems and the range of initial conditions leading 
to SNela explosions, together with other possible evolutionary effects (such as dependence 
on initial metallicity). The dependence on the time scale is illustrated in Figure 1 for the 
two main types of progenitor systems so far proposed, and for three different f)^{z). The 
systems labelled CLS(W) have evolutionary time scales of the order of a few Gyr whereas 
those labelled DD have time scales of the order of a few hundred million years only. The 
"flat" evolution of the efficiency for the DD systems for decreasing z is due to the fact 
that K/a(-2) closely follows f>*{z). Instead, £sNeia{,z) significantly increases towards lower 
z for the CLS(W) systems because the systems now exploding are a fixed fraction of the 
stars formed a few Gyr ago, when f>^{z) was higher. The fiattening of £sNeia{.z) at larger z 
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corresponds to a similar flattening of the f>^{z) considered. The absolute values of £sNeia{.z) 
reflect the abundance of progenitor systems (and thus the range of initial conditions leading 
to SNela), and they will be discussed below. 

The two types of systems whose SsNeia{z) is shown in Figure 1 have emerged in the last 
years as the main candidates to SNela progenitors. The DD systems are binaries made of a 
couple of CO WDs with a total mass exceeding the Chandrasekhar mass, and close enough 
to merge due to emission of gravitational waves in less than a Hubble time (Webbink 1984; 
Iben and Tutukov 1984). The CLS(W) systems consist of a CO WD plus a subgiant or 
red-giant companion that is overfilling its Roche lobe and transferring matter from its 
envelope to the WD {cataclysmic-like system or Algol-type binary that might be observed 
as a supersoft X-ray source). This scenario was originally proposed by Whelan & Iben 
(1973) and recently refined by Hachisu, Kato, & Nomoto (1996), who introduce the effects 
of a strong wind emitted by the accreting WD to stabilize mass transfer. That would allow 
the WD to reach the Chandrasekhar mass by steadily burning H into He and then He into 
CO (see also Nomoto, Iwamoto, & Kishimoto 1997). General discussions of the possible 
SNela progenitor systems can be found in Ruiz-Lapuente, Canal, & Burkert (1997). 

The SsNeia{z) values in Table 1 are based on p* = 2.0l!'lo;l8 x 10~^ /les Mq Mpc~^ yr^^ 
(Gronwall 1999) for the local SFR. For higher z, we give the average of the values calculated 
for each of the three different p*{z) referred to in Figure 1. The low-redshift efficiencies 
correspond to 1 SNela per ~ 9OOM0 going into star formation. As we see, SsNeiai^) is of 
the same order up to 2; ~ 0.55. Ongoing searches will soon yield the efficiency at z ~ 1. 
If we assume that all stars with initial masses M >,10 Mq produce SNell + SNelb/c 
(gravitational-collapse SN), then ^jj+ih/c{z) ~ 0.0054 x p*{z), for a Salpeter initial mass 
function (IMF) with x = 1.35, extending from 0.1 Mq up to 100 Mq. With the above 
values, that gives about 5 gravitational-collapse SN for every SNela in the local Universe 
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[see Tablel), in good agreement with the observations. 



We can relate £sNeia{z) to the fraction r] of stars in the 3 M© <jM <9 Mq initial mass 
range that should produce SNela if ^ia{z) closely followed p*(-2), since we would then have: 



J^^M)dM 
/o.fM$(M) dM 



£sNeIa{z) = r^ X ^00 ..^./..^ ... (2) 



where $(M) is the IMF. For the same IMF as above, we would have £sNeia{z) = 0.02304 x 77, 
which would mean r] ~ 0.06 for the {SsNeia) given below. Some 6% of the stars in the 
above mass range (making the approximation that all stars were born in binaries) should 
thus give rise to SNela. Given the strong selection that any initial population of binaries 
with primaries in that mass range undergoes in any scenario, at different evolutionary steps, 
until becoming a SNela candidate, the figure looks high if there is no help from the fact 
that most systems now exploding were formed at an epoch when p*(-2) was considerably 
higher than at the redshifts considered here, and it confirms the impression carried by 
Figure 1 that the DD systems, with their short evolutionary timescales adding to the strong 
evolutionary constraints, do fail to account for the bulk of observed SNela. We will see in 
the next Section that modeling of the Galactic DD population reinforces this conclusion. 
Comparing the measured EsNeia with that predicted from the DD scenario, at 2; ^ 0.55 
for instance, one sees that the prediction {EsNeia = 1-80 x 10~^ ^o^) is more than 2.95o" 
below the measurement {SsNeia = l-42to;42 x 10~^ ^o"*^)- The x^ test gives a probability 
P < 0.006 that the two efficiencies were the same. Even allowing for a sytematic error by 
a factor of 2 in the theoretical prediction of the rates and taking then the upper limit, the 
prediction would still be more than 2.5o" below the observational value (P < 0.012). If we 
refer to average values over the interval < z < 0.55, then we have, from observations, 
{SsNeia) = l-41^o:3? X lO"''^ M~^ whcrcas the DD prediction is 2.17 x 10"^ M-\ more than 
3.8CT below (P < 0.004). 



Precise enough measurements of ^jai^) at 2; >, 1 would in turn test the prediction that 
at those redshifts the decrease in metalhcity should inhibit the strong wind mechanism in 
the CLS(W) systems (Kobayashi et al. 1998). A decrease in £sNeia{z) should be observed, 
even if SNela do not completely disappear thanks to the early metal enrichment of elliptical 
galaxies. 

2.1. The local density of double degenerate binaries 

The accuracy of model predictions of the SNela rates expected from the merging of 
two CO WDs (the DD scenario) cannot just be tested against the actually measured rates, 
since the observed SNela can arise from a different evolutionary pathway. Only predicted 
rates much in excess of the measured ones would prove the modeling wrong. Since we find 
the efficency of SNela production in the DD scenario to be too low as compared with that 
observed, some extra test of the general binary evolution model is required. Such a test is 
provided by the measured space density of binary systems consisting of a pair of WDs (DD 
systems). In the DD scenario, a fraction of those systems (close CO WD pairs with a total 
mass exceeding the Chandrasekhar mass and close enough to merge in less than a Hubble 
time) would be the SNela progenitors. 

From existing surveys, no clear SNela progenitor candidate has been found among the 
detected DD systems, but aside from the role played by not too well controlled selection 
effects, the fact does not automatically validate model predictions that the space density of 
possible progenitors should be low enough for them to have escaped detection up to now. 
A better test is provided by the comparison of model predictions with the sample of DD 
systems already detected, even if no SNela candidate has yet been found. We have thus 
calculated, with the same Monte Carlo scenario code used to predict the SNela rates for 
different evolutionary scenarios (see Ruiz-Lapuente, Canal, & Burkert 1997, for instance). 
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the numbers of different types of DD binaries (He WD + He WD, CO WD + He WD, and 
CO WD + CO WD) that would form at times following an instantaneous burst of star 
formation. The formation of WD + MS pairs is also followed. That is later convolved with 
a SFR appropriate for the Solar neighbourhood (constant rate, ~ 1 Mq yr'^ for the whole 
Galactic disk), which gives us the number density of DD systems existing at the present 
time. Since WD pairs older than ~ 10^ yr should haved cooled below detection limits, the 
local SFR history is not really important when comparing the model with the observations. 
We find that the space density of WDs with ages below 10^ yr (wich roughly corresponds 
to My < 12.75) should be ~ 2.5 x 10"^ pc^^, which agrees well with the observational 
estimates of Liebert, Dahn, & Monet (1988) (~ 3 x 10"^ pc^^). 

In Table 2 (first line) we give the predicted fractions of close DD and close WD + MS 
pairs versus all WDs (single plus both wide and close pairs). We also give the fraction 
corresponding to the predicted SNela candidates. In the second line is shown the current 
birth rate of the three types of systems, and in the third one the numbers of systems with 
ages below 10^ yr. In the calculations we have adopted a value of the common envelope 
parameter acE = 1- Our results compare well with the observed samples (Saffer, Livio, & 
Yungelson 1998; Maxted & Marsh 1999). The predicted DD period and mass distributions 
do also agree (see Figure 2), especially if we take into account the observational biases 
against detection of systems with very low mass primaries (taken here to be the brightest, 
and thus the more recently formed members of the DD systems), and of binaries with either 
the longest or the shortest periods. Therefore, our modeling of the possible DD progenitors 
of SNela can be considered as being tested against all the currently available observational 
data, and our conclusion as to the low efficiency of the DD scenario in producing possible 
SNela progenitors, as compared with the measurements of the cosmic evolution of the 
SNela rates, appears robust. 
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3. Discussion and conclusions 

We have shown how, by using a cosmic approach, the topic of which stars end their 
evolution as SNela can be enlightened. About 6% of the binary stars in the range 3 — 9 M© 
end their lives as SNela. This fraction is unattainable from the merging of C+0 WD pairs, 
since only about 1% of binaries with masses in that range have a final total mass which 
is above the Chandrasekhar mass. The high production efficiency of SNela out of the 
star formation process found from 2; ~ to 2; ~ 0.5 and beyond suggests that the single 
degenerate scenario best explains Type la explosions. 

The above conclusion can be avoided if 

1. The star formation rate is understestimated by a factor of 6 even in the reconstructions 
from z ~ till z ~ 0.5 which give the highest estimates, while the rate of SNela is not 
underestimated. This situation seems unlikely as it would imply that the real SFR history 
is inconsistent with a dozen of empirical determination using very different methods. 

In addition, such extremely high SFRs would be in conflict with all we know about 
extragalactic background light. More especifically, at 2; = 0.55 there is just a discrepancy 
by a factor <, 1.5 among the SFRs in Figure 1. The bulk of DD systems merging at that z 
were formed at z ~ 0.62, where the range of possible SFRs is not any broader. 

2. We have overestimated by a factor of 6 the empirical rates of SNela from 2; ~ till 
z ~ 0.5. The possibility of having overestimated the SNela rate is unlikely since critical 
aspects such as obscuration by dust, or undetectability of SNela near galactic nuclei, should 
produce the opposite effect. 
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Tables 



Table 1: Type la supernova rates and production efficiencies along z 

Redshift tsnu ^la ^SNeia Probability 

(z) SNuh^s Mpc-V'his ^e^k (x' test) 

0} 0.2llg:f|| 2.2ttt 10-5 1-09+Jj^ 10-3 ** 

0.2 0.16±0.05 — — 

0.1^ O.Utom 1-7^1:? 10-5 1-24+^:8? 10-3 *** 

0.32^ < 0.32 (la) < 4.52 (la) 11.02 (2a) 10-^ < 2.27 - 5.54 10-^ *** 

< 6.2 (la) 15.0 (2a) 10-^ * 
0.38^ O.Sbtoil 4.8l|^ 10-^ 2.20ll;^} 10"^ 

6.9tti 10-5 * 
0.55^ 0.26± 0.08 4.53tl:i 10"^ l-42l[]:i 10^ 

6.74:tlll 10-5 * 

P < 0.004 
^Hamuy & Pinto (1999); ^Cappellaro et al. (1997); ^Hardin et al. (1999); ^Hamilton (1999); 
5Pain et al.(1996); ^Pain et al.(2000) *For the cosmology ^m =0.3 f]A=0.0; **SFR from 
Gronwall (1999); *** Averaging over the SFRs used in Figure 1; 
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Table 2: Predicted fractions of DDs and of WD + MS pairs 

Close DD CO + CO Pre-SNela CO + He He + He Close WD + MS 

1/15 1/114 1/419 1/24 1/58 1/6 

0.056 0.008 0.002 0.037 0.011 0.146 

5.6x10^ 8.5x10^ 2.2x10^ 3.7x10^ 1.1x10^ 1.5x10^ 
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Fig. 1. — Model predictions of SsNeia, the "efficiency" in producing SNela, that is the 
number of SN per unit mass of stars being formed, at a given z. The curves show its 
expected evolution for two SNela candidates systems with different timescales to explosion: 
merging of two CO WDs (dashed line, labelled DD), and cataclysmic-like binaries (or Algol- 
type systems), with the stabilizing effects of the accretion-induced wind included (solid line, 
labelled CLS(W), see text for further explanation). The results for three different star- 
formation histories are shown: Madau et al. (1998), Steidel et al. (1998), and Blain et al. 
(1998). The data points come from the SNela rate measurements of Pain et al. (1996, 2000), 
Hardin et al. (1999), and Hamuy & Pinto (1999). The results are independent from the 
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Fig. 2. — Top panel: Expected mass distribution of the brightest components of DD systems 
(sohd histogram). The dot-dashed histogram shows the distribution in the observed sample. 
Bottom panel: Expected period distribution of DD systems (sohd histogram). The dot- 
dashed histogram shows the distribution in the observed sample. The scarcity of low-mass 
WDs in the observed sample, as compared with the model prediction, can be explained by 
the increasing difficulty in detecting lower mass WDs, and the same applies to either very 
short or very long periods. 



